AO-AOO  M4 


UNCUAUIFXEO 


AD 

A0A3384 


NAVAL  OCEAN  SYSTCNS  CCNTCN  fAN  OXCtO  CALXF  F/«  17/2.1 
IXMULATXON  AND  MCASUHelCNT  OF  HXOM  FRCtUCNCY  XONOfFHCMC  CHANNC— CTCCU) 
AFR  77  D B SAILORS.  J R HILL 

NOSC/TR-111  NL 


NGSC 

• 

1 

1 

a 

1 

1 

1 

1 

eoNi 


9-77 


'V'i 


Technical  Report  111 


SIMULATION  AND  MEASUREMENT  OF 
HIGH  FREQUENCY  IONOSPHERIC 

CHANNELS 

Time-and  frequency-dispersive  properties  of  hf  sky 
wave  propagation  paths  can  be  modeled  for^use  in 
adapting  signal  format  and  signal  processing 
algorithms  for  communication  system  development. 

^/Sailors 

^ ! 4RyHill 

O ' i (Jj)  22  AprfI 


o ■' 


Interim  Beport.for  Period  July  1^6— Septe«toer  197^; 


Prepared  for 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
\ I PME  107-031 

! f ^ V I ; Washington.  DC  20360 


Approved  for  public  release;  distribution  is  unlimited  i - / 


/'  J/l 


NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CALIFORNIA  92152 


/ " / 


NOSC  / TR  111 


NAVAL  OCEAN  SYSTEMS  CENTER,  SAN  DIEGO.  CA  92152 


AN  ACTIVITY  OF  THE  NAVAL 

RR  GAVAZZI,  CAPT,  USN 

Commander 


MATERIAL  COMMAND 
HOWARD  L BLOOD,  PhD 

Technical  Director 


ADMINISTRATIVh  INFORMATION 

Work  was  performed  under  6272IN.  Task  XF2I22209I  (NOSC  B194B),  by  members 
o4  the  Propagation  Division.  This  interim  report  covers  work  from  July  through  September  1976 
and  was  approved  for  publication  22  April  1 977. 

The  authors  wish  to  express  their  indebtedness  to  ('C  Watterson  of  the  Institute  of 
Telecommunication  Sciences  (ITS),  who  sent  copies  of  ITS  reports  on  hf  channel  simulation 
and,  in  particular,  to  acknowledge  that  they  relied  heavily  on  the  ITS  report,  Hf  Channel- 
Simulator  Measurements  and  Performance  Analyses  on  the  USC-IO,  ACQ-6,  and  MX- 190  PSK 
Modems,  by  CC  Watterson  and  CM  Minister. 


Released  by 
JH  Richter,  Head 
Propagation  Division 


Under  authority  of 
PC  Fletcher,  Head 

Electromagnetic  Systems  Department 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whrn  Data  Enlarad; 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1 REPORT  number 

NOSC  Technical  Report  1 1 1 (TR  111)  / 

2 GOVT  ACCESSION  NO. 

3 RECIPIENT'S  CAT  ALOG  NUMBER 

4 title  SIMULATION  AND  MEASUREMENT  OF  HIGH 

FREQUENCY  IONOSPHERIC  CHANNELS.  Time-  and  frequency- 
dispersive  properties  of  hf  sky  wave  propagation  paths  can  be  modeled 

5 TYPE  OF  report  « PERIOD  COVERED 

Interim  — July  1976  — 

September  1976  ‘ 

for  use  in  adapting  signal  foniiat  and  signal  processing  algorithms  for 
communication  system  development. 

6 PERFORMING  ORG.  REPORT  NUMBER 

7 AuTmOR/^s; 

a CONTRACT  OR  grant  NUMBER''*) 

DB  Sailors,  JR  Hill 

9 PERFORMING  ORCANI  Z ATION  NAME  AND  ADDRESS 

Naval  Ocean  Systems  Center 
San  Diego,  California  921 52 

10.  program  ELEMENT.  PROJECT,  TASK 
AREA  a WORK  UNIT  NUMBERS 

62721N.XF2I222091 
(NOSC  B194B) 

1 1 CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Electronic  Systems  Command 

12.  REPORT  DATE 

22  April  1977 

(PME  107-031) 

13  NUMBER  OF  PAGES 

38 

1 U MONITORING  AGENCY  NAME  & ADDRESS^//  dltUrent  trom  ControtUng  OfftcB) 

15  SECURITY  CLASS,  (of  thia  report) 

Unclassified 

15«  DECLASSIFICATION  DOWNGRADING 

schedule 

16  distribution  statement  (of  this  Report) 

Approved  for  public  release;  distribution  is  unlimited. 

17.  distribution  statement  (of  the  ebetrsct  entered  In  Block  20.  H different  from  Report) 

IB  supplementary  notes 

' . I 

' ' 

19  KEY  WORDS  (Continue  on  reverse  side  If  neceeamry  end  Identityby  block  number. 

Modems  — high  frequency 

Simulation  - radio  reception 

High  frequencies  - distortion  , -A'-*' 

i'- 

lonograms 

Ionospheric  propagation  - refraction 
Wave  dispersion  — frequency  and  time 

20.  ABSTRACT  (Continue  on  reveree  aide  If  neceaaery  end  Identify  by  block  number) 

The  report  presents  literature  search  of  hf  channel  simulators,  describes  the  Watterson  model  and  its  limita- 
tions. and  discusses  distortion  of  ionosnhericallv  propaeated  signals  due  to  ionospheric  relMCtion. 

Significant  dispersion  occurs  if  the  half-bandwidth  exceeds  n l^y/  A r / A f,  where  At  / A l-(in  s/Hz)  is  the 
ionogram  slope  at  center  frequency.  For  100  kHz  bandwidth,  dispersion  is  significanCtf  At  / Af  exceeds 
3.14  X 10~4  ms/kHz,  its  measurement  requiring  62.8^^  resolution.  r . — 

Recommendations:  For  bandwidths  above  2.5  kHz  at  night  and  12  kHz  during  the  day,  add  frequency  dis- 
persion correction  filters  to  the  Watterson  simulator  tap  delay  lines.  To  achieve  required  time  delay  resolution, 
use  a complementary -sequence  or  modulated-pulse  oblique  sounder  system.  For  the  microprocessor  or  mini- 

DD  1473  EDITION  OF  1 NOV  65  |5  OBSOLETE 

S/N  0102  LF  014  6601 


UNCLASSIFIED 

security  CLASSIFICATION  OF  THIS  PACE  (Whmn  Data  KnlmrM) 


A 


UNCLASSIFIHD 

security  classification  of  this  PAGEflWien  Dmit  Enltrtd) 


Continuation  of  Block  20 

computer  controlled  receiver,  provide  an  atomic  time-frequency  standard  and  interfaces  for  terminal,  CRT  display, 
magnetic  tape  recorder,  and  auxiliary  devices.  When  measuring  modem  performance  capabilities  with  a channel 
simulator,  take  a “single  and  multiple  distortions”  approach.  Use  the  Jones  ray-tracing  program  to  obtain  the  fre- 
quency dispersion  relationships  needed  for  applying  corrective  filters  to  the  Watterson  model  ; relate  results  to 
channel-caused  system  distortions.  Use  a suggested  alternate  approach  if  too  many  delay  line  taps  would  be 
needed. 


UNCLASSIFIED 


security  classification  of  This  PAGEritTi»n  Dmtm  Enitrtd 


OBJECTIVE 


Establish  a channel  model  and  data  base  for  designing  sky  wave  modems.  In 
particular,  make  a literature  search  of  hf  channel  modeling  and  study  both  the  Watterson 
model  and  the  physics  of  the  channel.  Present  a plan  for  acquiring  available  information  to 
complete  any  necessary  modification  of  the  model. 


RESULTS 

1 . There  are  two  basic  sources  of  distortion  associated  with  hf  paths;  time 
dispersion  and  frequency  dispersion. 

2.  The  best  known  hf  simulator  is  that  developed  at  the  Institute  for  Telecommuni- 
cation Sciences  by  Watterson.  It  is  limited  to  system  bandwidths  less  than  2.5  kllz 
at  night  and  less  than  12  kHz  during  the  day. 

3.  An  important  source  of  distortion  of  ionospherically  propagated  signals  is  the 
dispersion  caused  by  ionospheric  refraction. 

4.  For  a 100  kHz  bandwidth,  dispersion  caused  by  ionospheric  refraction  is 
significant  if  Ar/Af  is  greater  than  or  equal  to  3.14  X 10'4  ms/kHz. 

5.  To  measure  dispersion  accurately  enough  to  provide  design  information  for 
dispersion  correction  filters,  62.8  ps  resolution  in  time  delay  is  necessary  lor  a 
100  kHz  bandwidth. 

6.  An  alternate  approach  is  presented  to  hf  channel  simulation  that  is  usable  when 
the  ionospheric  dispersion  is  so  large  that  it  would  call  for  an  excessive  number  ot 
taps  in  the  Watterson  model. 

RECOMMENDATIONS 

1 . When  using  the  Watterson  hf  channel  simulator,  install  frequency  dispersion 
correction  filters  on  the  tap  delay  lines  for  frequency  bandwidths  greater  than 
2.5  kHz  at  night  and  1 2 kHz  during  the  day. 

2.  To  achieve  the  needed  resolution  in  time  delay,  use  an  oblique  sounder  system 
employing  complementary  sequences  or  modulated  pulses. 

3.  Control  this  sounder  receiver  by  a microprocessor  or  a minicomputer. 

4.  Use  atomic  time  or  frequency  standards  for  timing  the  frequency  sweeps. 

5.  Use  a processor  that  has  interfaces  for  a terminal,  a CRT  display,  a magnetic  tape 
recorder,  and  an  auxiliary  device. 

6.  When  measuring  the  performance  capabilities  of  a modem  with  a channel 
simulator,  use  the  single-and-multiple-distortions  approach  discussed  in  this 
report. 

7.  Use  the  Jones  ionospheric  ray  tracing  program  to  determine  the  time  delay  versus 
the  frequency  and  range  dispersion  relationships  to  be  used  in  the  corrective  filters. 

8.  Write  a report  discussing  the  significance  of  the  various  system  distortions  that 
can  be  caused  by  the  channel.  Relate  these  results  to  those  obtainable  from  the 
Jones  ray  tracing  program. 


9.  When  the  spread  of  arrival  times  of  the  ionospherically  propagated  signals 
causes  an  excessive  number  of  taps  to  be  required  in  the  Watterson  model 
the  alternate  approach  to  ht  channel  simulation,  described  in  this  report. 
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INTRODUCTION 

I'wo  methods  are  eommonly  used  to  evaluate  the  performance  of  hf  radio 
communication  systems:  theoretical  analyses  and  experimental  measurements.  Both  arc- 
desirable  and  each  has  its  advantages  and  limitations. 

In  theoretical  analyses  of  system  performance,  various  hf  channel  models  have  been 
used.  Prediction  techniiiues  normally  used  are  a form  of  non-real-time  modeling  or  channel 
estimation  and  normally  do  not  take  account  of  the  more  transient  propagation  mechanisms. 

hxperimental  measurements  of  system  performance  have  been  made  over  actual 
ionospheric  links,  but  they  are  usually  difficult  to  obtain  and  to  apply.  Most  of  the 
difficulties  can  be  avoided  by  means  of  experimental  measurements  made  with  a stationary 
channel  simulator  in  the  laboratory.  Consequently,  there  has  been  a rapidly  growing  interest 
in  the  development  of  hf  channel  simulators.  (See  bibliography  entries  under  Watterson 
Model  and  all  but  the  last  two  under  Channel  Modeling.) 

An  hf  channel  simulator  provides  a means  for  a communications  engineer  to  optimize 
his  use  of  an  hf  propagation  path  that  has  dispersive,  time-varying  properties.  A simulator 
provides  access  to  the  values  of  parameters  of  an  appropriate  real-time  ionospheric  channel 
model  which  describes  the  path  behavior  adequately.  “Real  time”  in  this  context  implies 
that  the  parameter  values  can  be  updated  at  a rate  exceeding  the  maximum  significant  rates 
of  change  of  those  parameters.  Once  the  model  has  been  developed,  the  engineer  can  then  use 
it  for  optimally  adapting  the  signal  format  and  signal  processing  algorithms  for  communi- 
cation system  development. 

In  general,  there  are  two  basic  sources  of  signal  distortion  associated  with  hf  paths; 
time  dispersion  and  frequency  dispersion.  The  effects  of  these  mechanisms  can  be  represented 
by  the  “channel  scattering  function.”  This  is  a composite  function  that  shows  the  dispersive 
effects  of  the  path  on  an  ideal  impulse  in  the  time  domain  and  on  a single  cw  tone  in  the 
frequency  domain.  Distinct  propagation  modes  might  be  due  to  refraction  in  different 
ionospheric  layers,  multiple  refractions  in  individual  layers,  or  both.  Any  given  scattering 
function  is  a valid  description  of  the  path  behavior  only  for  a specified  rf  carrier  and  only 
while  the  “dispersion  surfaces,”  or  probability  distributions,  remain  sensibly  constant.  In 
the  longer  term,  time  variations  of  the  path  parameters  will  cause  changes  in  the  number  of 
modes,  the  mode  amplitudes  and  delays,  and  the  shape  of  the  distributions  associated  with 
each  of  the  modes. 

Although  an  unproven  channel  model  can  provide  useful  intormation,  results  must  be 
weighed  with  caution  since  they  may  not  be  typical  of  the  hf  medium.  More  useful  would  be 
a mathematically  tractable  channel  model  that  is  known  to  be  a valid  representation  of  the  hf 
medium  for  the  largest  practicable  percentage  of  typical  hf  links.  On  economic  grounds,  it  is 
important  for  the  channel  modeling  techniques  employed  to  be  no  more  complex  than 
necessary  and  to  be  tailored  to  provide  only  that  information  which  has  direct  relevance  to 
the  communicator.  As  an  example,  the  modeling  detail  needed  for  wide  band  systems  is 
greater  than  that  needed  for  narrow  band  systems  because  of  the  added  frequency  dispersion 
in  the  ionosphere.  For  narrow  band  systems  less  detail  would  be  necessary. 
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SYSTEM  DISTORTIONS 


In  any  radio  ct)innuinication  system,  the  sitmal  delivered  to  tlie  modem  re- 

eeiver  is  generally  a distorted  version  of  the  signal  that  the  modem  transmitter  should  ideally 
generate.*  The  soiirees  of  distortion  in  a digital  radio  com/mmications  system  are  as  follows: 

Transmitting  Equipment 

Modem  transmitter 
Rf  transmitter 

rhaiiiiel  (transmitting  to  receiving  antenna) 

Receiving  equipment 
Rf  receiver 
Modem  receiver 

The  types  of  transmitting  and  receiving  equipment  distortion  that  can  occur  are  as 

follows: 

Additive 

Thermal  noise 
Quantizing  noise 
Interference  ( internal ) 

Multiplicative 

Synchronization  (delay) 

Filtering  (frequency  selective) 

AGC  (time  .selective) 

Nonlinear 

The  types  of  distortion  that  a radio  channel  can  impose  on  a signal  are  as  follows: 

Additive 

Gaussian  noise 
Atmospheric  noise 
Manmade  noise 

Interference  (unwanted  signals) 

Multiplicative 

Multipath  time  scatter  (frequency  selective) 

Dispersion  time  scatter  (frequency  selective) 

Doppler  frequency  scatter  (time  selective) 

Fading  frequency  scatter  (time  selective) 


'office  of  Telecommunications  Report  75-56.  HF  Channel-Simulator  Measurements  and  Performance 
Analyses  on  the  USC-IO.  AC0-6,and  MX-190  PSK  Modems,  by  CC  Watterson  and  CM  Minister,  July  1975 
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All  tliese  types  ol' channel  distortion  are  present  on  ionospheric  communication  channels 
Dispersion  is  a variation  in  the  velocity  of  propagation  .with  frequency  over  any  path  m a 
channel.  It  causes  the  delay  characteristics  for  that  path  to  vary  with  fret)uency 

In  general,  the  small  distortions  introduced  hy  the  transmitting  equipment  are 
relatively  unimportant  and  can  he  ignored.  The  channel  distortions,  however,  can  he  suh- 
stantial  and  are  of  primary  importance.  While  the  distortions  introduced  hy  the  receiving 
equipment  may  not  he  larger  than  the  distortions  introduced  hy  the  transmitting  equipment, 
their  effects  on  the  system  performance  are  much  greater  when  multiplicative  distortion  is 
present  in  the  channel,  and  their  presence  cannot  he  ignored. 


WATTERSON  CHANNEL  SIMULATOR 

Watterson  et  al.,  of  the  Institute  for  Telecommunication  Sciences  (ITS),  undertook  a 
program  in  1965  to  develop  specialized  propagation  measuring  equipment,  to  obtain  measure- 
ments of  typical  ionospheric  channels,  and  to  analyze  those  measurements  to  determine  the 
validity  and  bandwidth  limitations  of  a proposed  channel  model.  Based  on  these  measure- 
ments, an  ITS  ionospheric  channel  simulator  was  built. This  simulator  is  presently 
capable  of  accurately  simulating  all  types  of  channel  distortion  listed  except  manmade  noise. 
They  found  that  their  channel  model  could  be  considered  valid  over  a bandwidth  that  is 
approximately  one-fourth  of  the  reciprocal  of  the  effective  time  spreads  on  the  ionospheric 
paths.  For  the  measured  channels,  these  bandwidths  were  2.5  kHz  for  the  nighttime  sample 
to  a maximum  of  1 2 kHz  for  the  daytime  samples.  They  feel  that  if  the  time  spreading  on 
each  measured  nighttime  path  were  caused  primarily  by  the  differential  propagation  time 
between  the  two  magnetic  components,  the  use  of  separate  taps  for  each  magnetoionic 
component  would  yield  a model  that  would  be  valid  over  a larger  bandwidth. 

A block  diagram  of  the  channel  simulator  as  it  is  used  in  non-space-diversity  measure- 
ments is  presented  in  figure  I . The  input  (transmitted)  signal  feeds  a delay  line  with  several 
taps,  one  for  each  propagation  mode  or  path,  numbered  1 .2.  ...  n,  where  n is  the  maximum 
number  of  modes  or  paths.  The  present  Watterson  simulator  shown  in  figure  1 has  four  taps. 
The  position  of  each  tap  can  he  independently  adjusted  to  a delay,  Tj.  This  adjustment 
sets  the  differential  path  delays  and  determines  the  multipath  distortions.  Each  tap  delivers 
a delayed  but  effectively  undistorted  version  of  the  transmitted  signal  to  a path  gain  multi- 
plier, where  the  signal  is  modulated  in  amplitude  and  phase  by  an  independent,  stationary, 
complex,  baseband  function  termed  the  path  gain  function,  (jj(t),  to  provide  the  frequency- 
scatter  distortion  (both  Doppler  and  fading)  on  that  path.  Each  path  gain  function 

Gj(t)  = (Gj  + Gj(t)|  exp  (j27rnjt),  ( I ) 


Vc  Watterson,  JR  Juroshek,  and  WD  Bensema,  Experimental  Confirmation  of  an  Hf  Channel  Model, 
IEEE  Trans  Commun  Technol,  v 18,  p 792-803,  December  1970 

^ESSA  Technical  Report  ERL  1 12-lTS  80,  Experimental  Verification  of  an  Ionospheric  Channel  .Model, 
by  CC  Watterson,  JR  Juroshek,  and  WD  Bensema,  July  1 969 

“^ESSA  Technical  Report  ERL  127-ITS  89,  Recommended  Specifications  for  Ionospheric  Channel  and 
Atmospheric  Noise  Simulators,  by  CC  Watterson  and  RM  Coon,  September  1969 
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Figure  1.  Clianiicl  simulator  for  noii-spaee-diversity  ineasurements. 
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where  Gj  is  a complex  nonfacling  specular  component  and  GjU ) is  a complex  fading  Gaussian- 
scatter  component  with  Rayleigh  amplitude  distribution  and  uniform  phase  distribution. 

Die  exponential  factor  introduces  a Doppler  frerjuetjcy  sliift,  n,-.  that  is  common  to  the 
specular  and  Gaussian-scatter  components. 

The  spectral  characteristics  of  each  patli  gain  function.  Gjtt).  are  described  by  the 
path  gain  spectrum.  Vj(n),  as  follows: 

V.(n)=z^  5(n-n:)+  7^'—  e'**' " . C) 

aj  V:7rajaj 


which  is  the  power  spectrum  of  Gjt  t ).  I'he  first  term  is  the  specular  component,  with  a 
Dirac-delta  function  spectrum  and  specular  attenuation  a-;  while  the  second  term  is  the 
Gaicssian-scatter  component,  which  has  a Gaussian-shaped  spectrum  with  a freciuency  spread 
( lading  rate)  o\  and  a Gaussian-scatter  attenuation  a'j.  1 he  specular  and  (iaussian-scatter 
attenuations  are  both  power  ratios:  the  transmitted  signal  power  divided  by  the  received  power 
in  that  component.  Fquation  (2)  is  illustrated  in  figure  2.  Five  parameters  together  completely 
specify  each  path  gain  function  and  path  gain  spectrum:  specular  component  attenuation, 
specular  component  phase.  Gaussian-scatter  component  attenuation,  Gaussian-scatter  com- 
ponent frec|uency  spread,  and  Doppler  sliift. 

Hie  delayed  and  modulated  signals  from  the  four  paths  in  figure  1 are  summed  with 
each  other  and  with  (iaussian  or  simulated  atmospheric  noise  to  form  the  output  (received) 
signal.  When  internally  supplied  cw  interference  or  other  externally  supplied  interference  is 
used,  it  is  introducetl  in  place  of  one  or  more  of  the  delayed  signals  so  that  independent 
Doppler  and  failing  can  be  introduced.  If  wide  band  externally  supplied  interference  isused. 
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Figure  2.  Path  gain  spectrum.*  i 

a second  delay  line  is  necessary  for  the  interference,  and  this  is  not  presently  available.  (Wide  i 

band  here  means  having  a bandwidth  larger  than  the  reciprocal  of  the  maximum  differential  f, 

delay  typical  of  hf.)  > 

In  hf  channel  simulation,  each  Gaussian-scatter  component,  Gj(t),  represents  in  f 

general  one  magnetoionic  component  when  the  two  magnetoionic  components  in  a high  or  | 

tow  ray  have  significantly  different  Doppler  shifts  or  delays.  When  the  Doppler  shifts  and  if 

delays  are  essentially  the  same,  a single  path  can  represent  both  components.  ’si 

The  channel  simulator  also  has  linear  dispersion  filters  that  can  be  inserted  in  an>  | 

path  between  the  delay-line  tap  and  the  path  gain  multiplier.  They  are  all-pass  filters  with  \ 

delay  characteristics  that  are  linear  in  frequency  and  with  selectable  amounts  of  dispersion  i: 

(delay  slope).  Since  the  filters  were  developed  for  vlf/lf  simulation,  they  do  not  necessarily 
represent  the  time-delay-frequency  slope  of  an  oblique  ionogram  and  are  not  usable  for  hi  ' 

simulation. 

Table  1 summarizes  the  numerical  values  of  the  channel  simulator  parameters  that 
are  available  in  the  ITS  Ionospheric  Channel  Simulator.  The  decibel  path  attenuations.  .\|.  jl 

are  defined  as  10  log  (aj).  The  attenuation  values  shown  are  accurate  over  the  range  ol  atten-  *! 

nations  normally  used  except  where  Aj  -*•  °°.  Here,  99  dB  settings  are  used,  and  accuracs  1 

is  of  no  practical  importance.  | 

fl, 


,ii 

•T 


,S 


TABLi:  1.  CMANNHL  SIMULATOR  SPKCIFICATIONS.' 


Parameter 

Available  Values 

Accuracy 

Bandwidth,  B 

3,  6,  or  12  kll/ 

— 

Number  of  Paths,  n 

4 maximum 

Delay,  7j 

0 to  10  ms  in  20-^s  steps 

0.5  ns 

Dispersion,  Dj 

0,  100,  220,  and  320/as/kll/ 

2 X 10-- 

Specular  At  ten,  Aj 

0 to  100  dB  continuous 

0.2  dB* 

Specular  Phase, 

continuous  (modulo  drr) 

0.01  radians 

Scatter  Atten,  Aj 

0 to  100  dB  continuous 

0.2  dB* 

Frequency  Shift,  v-^ 

Irrtl  Synthesizer:  0,  ±1 , ±2,  or 
+5  X 0.01, 0.1,  1,  10,  or  100  Hz 

1 X 10-4 

1 

Exi  Synthesizer;  -500  to  +500 
Hz  in  1 -mHz  steps 

1 X 10-b 

Frequency  Spread,  2oj 

Each  Path:  1,2,  or  5 XO.Ol, 
0.1,  1,  10,  or  100  Hz 

2X10-2 

One  Path.  1 .00,  1 .33,  1 .78, 
2.37, 3.16,4.22, 5.67, or 
7.50  XO.Ol, 0.1,  1,  10, or 
100  Hz 

2 X 10-2 

Fading  Spectrum,  Vj  (r) 

3-pole  approximation  of  Gaussian 

<0.03  dB** 

Signal  to  Noise,  Eij/N'j 

-100  to  +40  dB  continuous. 

Gaussian 

0.2  dB* 

Atm 

0.5  dB* 

Atmospheric  Noise 

hf:  7.2  dBm  2.7  kHz 

— 

Impulsiveness,  Vj 

If:  10.8  dB  in  2.7  kHz 

— 

If:  9.0  dB  in  2.7  kHz 

Signal  to  Interference,  S/I 

-1 00  to  +90  dB  continuous 

0.2  dB* 

Space  Diversity 
Internal  Distortions 

Single  or  independent  dual 

— 

Additive 

<-70  dB  relative  to  signal  + noise 

— 

Nonlinear 

<-70  dB  relative  to  signal  + noise 

*Over  range  of  significance 

* ‘Theoretical 

CHANNEL  CHARACTERIZATION 


i 

\ 

» 

i 


In  analyzing  the  performance  capabilities  of  modems  suffering  multiplicative 
channel  distortions,  it  is  useful  to  have  theoretical  statistical  descriptions  of  the  channel 
frequency-scatter  and  time-scatter  characteristics  in  addition  to'the  channel  simulator 
parameters.  Two  statistical  descriptions  for  stationary  time-varying  channels  given  by 
Bello^  and  others  are  the  channel  correlation  function  and  the  channel  scatter  function, 
which  can  be  obtained  from  the  channel  time-varying  frequency  response.  From  figure  1 
and  equation  ( 1 ),  it  can  be  seen  that  the  complex  channel  time-varying  frequency  response 
for  the  ITS  channel  simulator  model 

n 

H(f.t)=  2 Gj(t)  exp  t-jZTTTjf)  (3) 

i=l 


and  the  channel  correlation  function 


lim  ' 

R(Af,At)=."!!^  — / H*  (f,t)H(f+Af,  t + At)  dt  , 
M “M  -t 


(4) 


where  the  asterisk  (*)  denotes  a complex  conjugate.  The  channel  scatter  function  then  is 
the  double  Fourier  transform  of  R (Af.At): 


-oo  -oo 

S(r,.)=y  J 

—oo  _oo 


R(Af,At)  exp  (j27rrAf-  InrAt)  dAfdAt. 


(5) 


For  the  ITS  channel  simulator  model,  it  can  be  shown  that  the  channel  scatter  function 


n 

S(T,r')=  2 5(r  - Tj)  Vj!!-)  . (6) 

i=l 

where  Vj(t')  is  the  path  gain  spectrum  specified  by  equation  (2).  While  S(T,n)  is  dimension- 
less. it  is  a power-ratio  density  function;  it  specifies  the  power  in  the  received  signal  per 
unit  delay  and  per  unit  frequency  offset,  divided  by  the  transmitted  power,  as  a function  of 
delay,  r,  and  frequency  offset,  n.  It  specifies  how  the  signal  is  scattered  or  distributed  in  time 
by  the  multipath  delays  and  in  frequency  by  the  Doppler  shifts  and  fading.  It  is  not  a des- 
cription of  a signal,  but  rather  what  the  channel  does  to  a signal. 

Additional  statistical  descriptions  are  also  useful  (Watterson,  et  al)'^.  The  channel 
frequency-scatter  function 

yjv)  = f Str.i’jdT 

'Loo 


n 

i=] 

^PA  Bello,  Characterization  of  Randomly  Time-Variant  Linear  Channels,  IEEE  Trans  Commun  Syst  1 1 , 
p 360-393,  December  1963 
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riie  diannd  fraiuency-scattcr  function,  like  the  path  gain  spectrum  of  each  path,  can  be 
described  in  terms  of  a Doppler  shift  and  frequency  spread.  To  do  this,  it  is  necessary  to 
normalize  equation  (7)  so  that  its  integral  is  unity.  Since  S(t,j^)  is  a power-ratfo  density 
function,  it  follows  that  the  channel  attenuation 


1 


Str.r’)  drdt^ 


(Si 


and  the  normalized  channel  frequency-scatter  function  isa^.V^.(r').  The  channel  Doppler 
shift  then  is  the  first  moment  of  the  normalized  channel  frequency-scatter  function; 


/ 


la^.V^.t^ll 


di' 


(^)) 


where  equation  (2)  and  equation  (7)  were  used  to  obtain  the  last  form.  I he  channel  tre- 
quency  spread,  2o^.,  is  defined  as  twice  the  square  root  ot  the  second  central  moment  ot 
the  normalized  channel  frequency-scatter  function: 


1/2 


n 

f \ 

ur  i>r  + a- 

^c  2 

. i=l 

11  1 

I/: 


(10) 


These  definitions  of  channel  Doppler  shift  and  channel  freciuency  spread  are  consistent  with 
the  corresponding  path  definitions;  ie  for  a single  path  channel,  equations  (‘>|  and  (10) 
yield  the  path  Doppler  shift  and  path  frequency-spread  values  of  equations  ( 1 ) and  ( 2 1. 

The  time-scatter  characteristics  of  the  channel  can  be  described  in  a similar  way. 

The  channel  time-scatter  function 


(11) 


i 


8 


riic  channel  delay,  is  ihe  first  moment  of  the  normalized  channel  time-scatter  function: 


and  the  channel  time  spread,  2p^..  is  twice  the  square  root  of  the  second  central  moment  of 
the  normalized  channel  time-scatter  function; 


If  time  in  the  preceding  equations,  ( 1 ) through  ( 1 3),  is  measured  relative  to  the  transmitter 
timing,  all  path  delays,  Tj,  are  positive,  and  the  channel  time  delay  is  usually  much  larger 
than  the  channel  time  spread.  Since  a constant  time  delay  can  be  subtracted  from  the  channel 
(ie,  the  same  constant  subtracted  from  all  path  delays)  without  changing  the  distortion 
present  on  the  signal,  it  is  convenient  in  analyzing  modem-receiver  performances  to  measure 
time  such  that  the  channel  delay  in  equation  (12)  is  zero.  When  this  is  done,  the  channel 
scatter  function  in  equation  (6)  for  a channel  with  two  Gaussian-scatter  paths  (no  specular 
components)  would  appear  in  general  as  shown  in  figure  3.  The  two  paths  might  represent 
one-  and  two-hop  propagation  modes  via  the  F layer  of  the  ionosphere.  Because  the  channel 
time  delay  was  set  to  zero,  the  first  path  has  a negative  delay  T^  and  the  second  path  has  a 
positive  delay  T2- 
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Figure  3.  Scalier  lunclioii  tor  a iwo-palh  channel. 

DISPERSION  DUE  TO  IONOSPHERIC  REFRACTION 


A source  of  distortion  of  ionospherically  propagated  signals  is  the  dispersion  as  a 
consequence  of  ionospheric  refraction.  This  dispersion  manifests  itself  hy  altering  the 
phases  of  the  various  frequency  components  of  a signal  so  that  they  are  not  coherently 
combined  by  the  correlator.^’ 

The  phase  distortion  imposed  on  a signal  by  ionospheric  dispersion  may  be  repre- 
sented by  a transfer  function  of  the  form  * where  the  argument  0(f)  represents  thc 

variation  with  frequency  of  the  phase  retardation  over  the  propagation  path.  If  0(f)  is  a 
linear  function  of  f.  there  is  no  dispersion,  and  the  transfer  function  then  represents  only 
the  propagation  delay  over  the  path.  Such  a group  delay,  r.  may  be  represented  as 

d0(  f ) , , 

T= , a constant.  (14) 

df 

However,  if  the  variation  of  0 with  frequency  is  nonlinear,  so  that 
d0(f) 

T = , not  a constant , ( 1 e ) 

df 


then  dispersion,  with  its  conseciuent  distortion,  is  present. 

^Klectrophysics  Laboratories  Project  Report  236,  Dispersion  Correction  Study,  by  RN  DeWitt.  JR  Earhart. 
and  JD  Kolesar,  August  1973 

^Electrophysics  Laboratories  Project  Report  243.  Dispersion  Correction  Study,  by  RN  DeWitt,  JL  Wheeler, 
and  FJ  WOods  Jr.  March  1974 

*^Eleclrophysics  Laboratories  Project  Report  207,  Estimates  of  Distortion  of  a Wideband  Signal  Propagated 
via  the  Ionosphere,  by  R.\  DeWitt.  JR  Earhart,  and  JD  Kole.sar,  .September  1972 
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In  the  development  of  the  Watterson  tap  delay  simulation  mode,  it  was  assumed 
that  the  time  delay  was  stationary  with  frequeney.  Watterson  himself  found  that  his  model 
was  limited  to  a bandwidth  of  I 2 kHz  for  daytime  samples  and  2.5  kHz  for  nighttime  sam- 
ples. If  this  model  is  to  be  used  for  larger  bandwidths,  it  must  incorporate  some  means  for 
tolerating  or  adapting  to  the  existing  ionospheric  phase  shift  and  delay. 

An  examination  of  ionospheric  sounding  seems  to  indicate  that  the  group  retardation 
over  the  path  varies  smoothly  with  the  portion  of  the  hf  spectrum  where  any  given  propaga- 
tion mode  exists.  Tor  a fixed  path,  the  group  delay  distortion  may  be  different  for  each 
mode.  If  such  is  the  case,  corrective  filters  could  be  added  to  each  tap  or  precorrection  of 
the  signal  could  be  feasible.  The  amount  of  dispersion  on  a signal  is  a function  of  the  partic- 
ular ray  paths  between  the  terminals  of  the  communication  link.  Consequently,  the  disper- 
sion will  vary  as  the  path  end  points  vary,  and  a correction  made  for  a signal  between  the 
fixed  tenninal  and  a remote  point  is  not  appropriate  to  a path  between  the  fixed  terminal 
and  another  point  sufficiently  displaced  from  the  first. 

In  propagating  via  the  ionosphere,  if  the  signal  experiences  a frequency  dependent 
phase  delay  0(f),  we  may  approximate  0(f)  by  a Taylor  series  expansion  about  the  band 
center,  f^,.  Thus 

0(f)  = 0(i;,)  + 0;,(fo)(f-to>+ '/-^o**o)<‘'-‘o)“  + 


We  ignore  the  first  two  terms,  since  the  Watterson  model  includes  these  terms  in  equation 
(3).  Assuming  that  the  third  term  well  represents  the  remaining  variation  of  phase,  the 
variation  may  be  regarded  as  seriously  degrading  to  the  signal  if,  over  the  interval  extending 
from  the  band  center  to  the  band  edge  - a difference  of  b/2  in  frequency  the  phase  con- 

tribution of  the  third  term  varies  through  a half-cycle.  That  is,  if 


1/2 


(b/2)-  = 


7t/2 


(17) 


Solving  for  0',  (f^,),  we  find  that  this  condition  occurs  if 


A,**  ! f I — 


47t 


(18) 


The  dispersion  may  be  regarded  as  significant  if  b/2  > 


, where  Ar/Af  in 


s/Mz  is  the  slope  on  the  ionogram  at  f^,  the  center  frequency.  If  the  spectral  components 
are  summed  without  correction  under  such  circumstances,  the  components  of  the  signal  at 
band  center  will  tend  to  be  canceled  by  the  frequency  components  of  the  signal  at  band 
edge.  Intermediate  components  will  be  only  partially  canceled;  thus  the  condition  does  not 
imply  a strong  degradation  of  the  signal.  Rather  it  represents  the  point  at  which  dispersion 


begins  to  become  serious.  F'igure  4 is  a plot  of  b/2  in  kHz  equal  to  • where  At  Af 

is  in  ms/kllz.  If  the  dispersion  has  been  corrected,  figure  4 is  a measure  of  the  resolution 
necessary  to  obtain  adequate  dispersion  data  for  simulation  purposes.  At  b/2  of  50  kHz.  for 
instance,  the  resolution  necessary  would  be  62.8  /as  in  time  delay. 


J 


F'igure  4.  Ionospheric  bandwidth  dispersion  limit. 

Figure  5 is  a backscatter  ionogram  I'roni  the  Chapel  Bell  Radar  at  Whitehouse.  Vir- 
gitiia.^  The  slope  at  the  t'rei|uency  of  measitremetU  was  5 X 10“^  ms/kHz,  This  corresponds 
to  a distortion  half-band,  using  the  above  criteria,  of  260  kHz.  Straight  lines  are  drawti 
tangent  to  the  oblique  ionogram  traces  at  the  frequency  of  measurement,  for  two  targets. 

A line  parallel  to  the  upper  trace  is  drawn  so  as  to  intersect  the  lower,  and  the  difference  in 
slopes  of  the  lines  is  determined  from  a measurement  of  their  separation  at  a distance  from 
their  intersection  corresponding  to  3 MHz.  This  difference  divided  by  the  time  separation 
of  the  two  ionogram  traces  at  the  frequency  of  measurement  gives  the  range  dependence  of 
dispersion. 

At  the  MUF-  of  any  of  the  received  modes  dr  'df  = <».  If  the  upper  edge  of  the  band- 
width b is  at  the  MUF,  a worst  case  for  dispersion  is  established.  The  boundary  condition 
defined  by  the  MUF  can  l>c  used  near  the  MUF  to  define  a power  series  fit  to  the  ionogram 

n+l 

j=l 
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Figure  5.  BackscuUer  io/iogra/ii  1‘rom  Wliilehouse  VA.  scaled  to  deter- 
mine  the  range  dependence  of  dispersion  ( 10  14/71.  I 252Z).^ 
Note:  A backscatter  ionogram  shows  reflections  for  several 
ranges;  whereas  an  oblique  incidence  ionogram  is  for  only  one 
range,  liach  range  has  a different  slope. 


BEARING  290 


FREQUENCY,  MHz 


where 

T = delay  time  (ms) 
fj,  = carrier  freiiiieney  (MU/) 
f = operating  fretiiieney  (.Mil/.) 
n = order  t)f  polynomial 
Aj  = power  series  eoelfieients. 

1 or  the  power  series  model,  the  phase  distortion 

n+2 
A<>,  = 2 

J=3 

It  IS  awkwanl  to  work  with  a power  series  expansion  in  this  respeet  because  an  a 
ment  ol  the  eoeffieient  of  a higher  term  to  improve  the  fit  will  affect  the  correctness  ol  any 
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prior  adjustiiKMU  of  lower  order  eoelTieieiits.  To  avoid  this  interaction,  it  is  preferable  to  j 

employ  a set  of  orthogonal  functions,  the  adjustments  for  which  are  independent  of  one  j 

another.^  Orthogonalization  of  simple  power  laws  over  the  bandwidth  of  the  signals  lends  || 

to  phase-correcting  functions  which  are  Legendre  polynomials/’ 

The  amount  of  disix-rsion  appearing  on  the  oblique  ionogram  varies  with  time  of 
day,  path  range,  season,  atid  mode  as  well  as  witii  environmental  conditions.  Usually  the 
slope  will  be  larger  at  night  than  during  the  day  and  at  short  ranges. 

EXTENSION  OF  SIMULATION  TECHNIQUE  TO  WIDE  BANDWIDTHS 

The  time  delay  for  each  mode  or  propagating  path  is  a measurable  quantity  over  the 
frequency  bandwidth.  If  we  assume  that  the  time  delays  for  the  i**^  mode  are  known  as  a 
function  of  frequency  and  range,  then  the  time  delay  can  be  approximated  by  a polynomial 
function  of  frequency  and  range: 

n+l  j-H 

^i"o-Ko>=  S S A^jtf-f^f^-ltR-Rofi-'.  CD 

11=1  j=l  I 

where  f,,  and  R,.  are  reference  frecpiency  and  range.  The  coefficients  can  be  determined 

1)  t)  Q 

from  measurements  or  trom  ionospheric  ray  tracing.'' 

Having  determined  r(f,t).  then  a new  Gj(t,f)  is  found  by  using  the  relationship 

Gj(t,f)  = gj(t)+gj(f).  (22) 

I quation  ( 22)  includes  gj(t),  which  is  the  Gj(t)  of  the  Watterson  model. 

The  function  of  gj(f')  can  be  applied  by  adding  a corrective  filter  to  each  tap  which 
would  represent  the  dispersion  spectra,  g|(f),  relative  to  fy,  the  bandwidth  center.  Each  fy 
would  have  a different  filter.  The  number  of  taps  would  depend  on  the  amount  of  correc- 
tion desired  for  each  filter  and  would  be  equal  to  the  bAT  product  of  the  ionospherically 
distorted  signal,  where  b is  the  bandwidth  and  AT  is  the  range  of  times  over  which  the 
various  modes  are  spread. 

Typically,  the  spread  of  arrival  times  might  be  a few  milliseconds,  so  for  a 100  kHz 
signal  a few  hundred  taps  would  be  required.  One  could  also  model  the  various  modes 
individually,  with  fewer  taps.  For  example,  if  the  F-scan  of  the  path  indicates  a mode  has 
a range  of  delay  of  AT  over  the  bandwidth  of  interest,  W,  one  may  be  able  to  make  a satis- 
factory signal  mode  model  with  WAT  taps.  However,  it  may  prove  to  be  awkward  to  orches-  I 

trate  several  separate  modes  in  this  fashion.  n 


^^Office  of  Teloeommunieations  Report  75-76,  A Versatile  Three-Dimensional  Ray  Tracing  Computer  Pro- 
gram for  Radio  Waves  in  the  Ionosphere,  by  RM  Jones  and  JJ  Stephenson,  October  1975.  Tllis  is  a versa- 
tile program  for  tracing  rays  through  an  anisotropic  medium  whose  index  of  refraction  varies  continuouslv 
in  three  dimensions.  The  program  can  represent  the  refractive  index  by  either  the  Appleton-Hartree  or  the 
Sen-Wyller  formula.  The  program  has  nine  ionospheric  models  for  election  density,  eight  perturbations 
(irregularities)  to  the  electron  density  models,  four  models  of  the  earth's  magnetic  field,  and  four  models 
of  the  electron  collision  frequency. 


I 
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I'or  a real-time  model,  tapped  delay  lines  as  discussed  above  might  be  well  suited  to 
the  task.  The  necessary  tap  weights  could  be  computed  fast  enough  in  a small  processor  to 
keep  them  suitably  updated  for  modeling  gradual  ionospheric  variations,  while  the  iono- 
spheric “filter”  would  take  place  in  the  time  domain  and  could  take  place  at  a much  higher 
rate.  A similar  process  could  be  used  to  correct  distortion  in  a working  system. 

AN  ALTtRNATE  APPROACH  TO  HF  CHANNEL  SIMULATION 

RN  DeWitt'*^  has  expressed  some  reservations  about  the  Watterson  channel  simu- 
lator. These  mostly  concern  his  use  of  tap  weights  that  are  “complex  (bivariate)  Gaussian 
stationary  ergodic  random  processes”  with  zero  mean.  IX-Witt  feels  that  this  kind  of  model- 
ing seems  to  be  more  of  a carryover  from  early  low  resolution  studies.  That  is,  ionospheric 
signals  have  always  been  observed  to  undergo  a great  deal  of  nuctuation  as  a result  of  unre- 
solved multiple  propagation  modes  that  interfere  with  one  another. 

However,  there  is  some  ciuestion  as  to  whether  these  same  fluctuations  would  appear 
on  a single  resolved  mode.  Tor  the  modes  to  be  entirely  resolved,  a rather  wide  bandwidth 
would  be  required,  particularly  in  the  case  of  magnetoionic  splitting.  But  once  the  modes 
are  resolved,  what  mechanisms  exist  to  produce  fluctuations?  Of  course,  there  are  some.  For 
example,  there  are  absorption  effects,  such  as  SIDs,  that  will  reduce  signal  strength,  and  TlDs. 
which  may  either  increase  or  decrease  signal  strength  of  a given  mode  through  focusing. 
Diffraction  by  sporadic  E irregularities  could  constitute  a third  mechanism.  These  phenom- 
ena are  not  always  present,  however,  and  even  when  they  are,  it  is  unlikely  that  they  are  well 
modeled  by  a bivariate  Gaussian  zero  mean  process.  Diffraction  from  irregularities  is  the  one 
that  could  sometimes  produce  a Gaussian  process. 

An  alternate  approach  to  hf  channel  simulation  is  presented  here.  This  approach 
provides  the  transfer  function  for  a single  propagation  mode;  the  effects  of  multimoding 
would  have  to  be  handled  by  a superposition  of  such  transfer  functions.  A subset  of  a 
receiver  system  has  been  described.  The  character  decoding  can  be  done  many  other  ways 
but  has  been  included  here  to  illustrate  the  use  of  the  ionospheric  model  applied  to  both  a 
random  waveform  and  a pure  character  waveform. 

A waveform,  W(t).  can  be  written  as  a modulation  function.  M(t),  and  a carrier  wave, 
cos  ( 27rf^.t  + <p^.): 

W( t ) = M( t ) cos  ( 27rf^.t  + 0^.)  . ( 23 ) 

The  frequency  independent  portion  of  the  function  Gj(t.f)  can  be  included  in  the  model 
when  each  sample  W(t)  is  calculated.  Thus 

Wj  (t  +gj(t))  = M(t  +gj(t))  cos(27Tf^,(t  -i-gj(t))  + 0^.)  . (24) 


'^International  Telephone  and  Telegraph  Avionics  Division  Electro-Physics  Laboratories  Letter  serial  no. 
RD-H.186,  by  RN  DeWitt.  to  Naval  Electronics  Laboratory  Center  (DB  Sailors).  Subicet : Review  of 
Watterson ’s  IEEE  paper,  ,)0  September  D)7(r 
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wiiATC  gj(t)  is  a ranclt)m  function  having  Die  desired  characteristics,  eg  (Jaiissiaii,  nonstatiunary, 
etc.  The  fjtt  ) portion  ot  Cij  is  folded  into  the  model  alter  taking  the  last  l oiirier  transtorm 
( FT  T)  of  Wjtt ).  The  spectrum  of  the  final  waveform  W'(t ) 


nTlW'itM  = n-TlWtt  +g| 


I f 

-2ni  T:(n)di^ 

>o  ' 


(25) 


Niite  that  if  r is  represented  by  a polynomial  in  f as  in  (21  ).  then  the  exponential  in  ( 25) 
and  (3)  is  a phase  shift 


^■(f)  = 0j(f)  - 27rf(aj  + bjf/2  + Cjf-.'3  + . . . ) 


(2b) 


DECODING  THE  CHARACTER  FROM  THE  SPECTRUM 

Correlations  can  lx-  calculated  very  rapidly  in  the  freciuency  domain.  We  are  nearly 
through  at  this  point,  since  we  already  have  the  I E F ol  W'.  By  using  precalculated  stored 
s|x-ctra  of  the  character  set  we  have  immediately  the  correlation  function  S(^  of  the 
character; 


rrris^(t)i  = i i ric'tt)!  x iTT|w'(t)i  . 


(27) 


The  function  S^(t)  will  peak  above  the  threshold  at  a value  of  t corresponding  to  the  posi- 
tion of  the  character  in  the  original  wave. 


SPECTRA  OF  THE  CHARACTER  SET 

The  spectrum  of  each  character,  M'T(C'(t)l , is  calculated  by  using  the  same  process 
as  that  of  W"(t).  That  is,  each  of  the  modes  included  in  the  calculation  of  W'  uses  the  time 
delays  t.  The  only  difference  between  C'(t)  and  W'(t)  is  that  gj(t)  is  set  to  zero  for  calculat- 
ing C'(t).  This  is  because  the  mean  time  delay  (nonrandom  part)  is  accounted  for  in  t. 

Tims  the  F'FT(C')  can  be  precalculated  when  Tj(f)  is  measured. 

SAMPLING  CONSIDERATION 

Assume  that  a signal  with  bandwidth.  B,  is  processed  in  a linear  filter.  No  signals 
higlier  in  fretiuency  than  f^.  + B nor  lower  than  f^.  - B will  be  produced.  If  B = 100  kHz,  the 
sampling  rate  must  be  at  least  once  every  5 seconds  if  the  Nyqiiist  frequency  , fjq,  is  to 
equal  B.  1024  samples  will  cover  about  5 ms,  which  should  be  enough  for  forward  and 
inverse  transforms  to  be  calculated  in  real  time.  The  resulting  frequency  resolution  is  200 
llz. 
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ALIASING 


The  earriei  frequency  can  be  “shifted"  to  zero  frequency  provided  that  the  Nyquist 
frequency  (sampling  interval)  is  related  to  the  carrier  Irequency  as  lollows: 

fj^.  = 4mf[s^.  where  m is  an  integer.  < 2H) 

The  transform  will  be  shifted  in  freriuency  so  that  I'  = t - 1^.. 

N ^ ' c " * N 

0-f,  <2^) 

+*  N '"c  + ' N 

The  requirement  here  is  that  the  incoming  waveform  is  band-limited  to  plus  and  minus  B. 


VERIFICATION  OF  THE  CHANNEL  MODEL 

Preferably,  both  for  theoretical  analyses  and  as  the  basis  for  simulator  designs,  a 
desirable  channel  model  would  be  one  that  is  known  to  be  a valid  representation  of  the  hf 
medium  for  the  largest  practicable  percentage  of  typical  hf  links.  Some  of  the  detailed 
characteristics  of  the  ionosphere  needed  to  specify  such  a model  may  not  have  been 
measured. 

In  1 969  Watterson  reported  experimental  verification  of  his  ionospheric  channel 
model. A hybrid  time-and-frequency  domain  method  was  used  to  obtain  propagation 
measurements.  In  the  transmitter,  an  accurate  cesium-beam  frequency  standard  was  used 
to  generate  a carrier  that  was  amplitude  modulated  to  obtain  coherent  pulses  with  near- 
Gaussian  envelopes  of  50  /as  duration  at  a ZOO-Hz.  rate.  The  same  reference  was  used  to 
generate  1 1 cw  signals  (tones)  of  equal  amplitude  at  frequencies  that  coincided  with  1 1 
spectral  components  of  the  pulses.  The  transmitter  output  (both  the  power  spectrum  and 
the  time  function)  is  illustrated  in  figure  6.  Tlie  1 1 tones  were  symmetrically  spaced  in  a 
9.6  kHz  band  about  the  carrier  frequency,  f^,.  The  low  level  composite  signal  was  amplified 
linearly  to  5 kW.  In  the  receiver,  each  of  the  1 1 tones  (with  delays  and  modulation  imposed 
by  the  ionospheric  channel)  was  separated  and  heterodyned  to  a second  IF  of  80  Hz.  Ttie 
received  signal  was  also  pas.scd  through  an  1 1 -notch  filter  that  removed  the  1 1 tones,  and 
the  resulting  pulses  were  envelope  detected  and  compared  against  a 200  Hz.  reference  to 
obtain  measurements  of  t\.  The  80  Hz.  IF  tones  were  FM  recorded  on  an  analog  tape- 
recorder  with  an  80  Hz.  reference  signal.  During  playback,  the  reference  signals  were  used 
to  detect  each  of  the  1 1 IF  tones  to  obtain  the  real  and  imaginary  components  of  the  time- 
varying  frequency  responses  of  the  measured  channel,  Hpf,(f|^,t).  at  the  1 1 values  of  f|^. 

Propagation  measurements  were  made  over  a 1 294  km  link  between  Long  Branch. 
Illinois  and  Boulder,  Colorado.  Nighttime  measurements  were  made  at  5.86  MHz  and  day- 
time measurements  at  9.259  MHz,  both  far  enough  below  the  predicted  MUF  that  both  one- 
and  two-hop  modes  probably  would  be  seen.  About  7 hours  of  recorded  nighttime  measure- 
ments and  10  hours  of  daytime  measurements  were  obtained  during  three  periods  in  Novem- 
ber 1 967.  Times  were  avoided  near  sunrise  and  sunset,  when  propagation  conditions  were 
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clKinginj;  rapidly.  Samples  were  selected  that  appeared  typical  and  reasonably  stationary. 
Three  of  these  samples  of  10  to  13  minute  duration  from  day  and  night  measurements  were 
analyzed  to  verify  the  model.  The  limiting  bandwidths  were  1 2 and  8 kHz  on  the  two  day- 
time samples,  2.5  kHz  on  the  nighttime  sample.  Since  the  1 1 tones  were  symmetrically 
spaced  in  a 9.6  kHz  band  about  the  carrier  frequency,  it  is  difficult  to  interpret  the  eflects 
of  bandwidth  beyond  9.6  kHz.  Nor  are  three  samples  sufficient  to  state  that  they  are  a 
valid  representation  of  the  hf  medium.  The  sunrise/sunset  data  would  have  provided  valua- 
ble data  on  dispersion  if  the  dispersion  alone  had  been  analyzed.  Since  a sounder  was  not 
run  coincidentally  with  the  experiment,  it  was  not  known  how  far  the  observation  Irequen- 
cies  were  from  the  operational  MUF. 


CHANNEL  MEASUREMENTS  AND  THEIR  PRACTICAL  IMPLEMENTATION 

Attention  will  now  be  turned  to  three  basic  classes  of  channel  measurements  which 
require  the  use  of  special  purpose  equipment. 


CLASS  1 CHANNEL  MEASUREMENTS 

A channel  stimulus,  x(t),  is  transmitted  repetitively  over  the  hi  channel  from  a base 
station.  Apart  from  noninterference  with  other  spectrum  users,  there  is  no  constraint  on  the 
form  of  x(t);  it  could  be  a relatively  complex  coded  signal.  Specialized  processing  equip- 
ment at  a remote  station  analyzes  the  received  version,  y(t),  of  the  channel  stimulus  and 
formulates  an  appropriate  path  model.  In  the  modeling  process,  channel  estimation  data 
from  other  sources  can  be  incorporated  legitimately  if  they  are  appropriate. 

The  best  known  implementation  of  Class  I channel  measurements  is  oblique  incidence 
ionospheric  sounding,  of  which  several  forms  are  available.  The  simplest  type  is  a pulse 
sounder  employing  time  and  frequency  synchronized  transmission  and  reception.  A high- 
power  sounding  transmitter  radiates  short  pulses  at  given  times  on  preset  channels  that  arc- 
spaced  in  frequency  over  the  entire  hf  band  or  a part  of  it.  Tlie  transmission  schedule  can 
be  varied  as  required  by  a system  “program”  and  is  commonly  controlled  via  a clock  which 
can  be  locked  to  an  external  time  standard  transmission  such  as  WWV.  At  the  remote  station, 
the  sounder  receiver  is  synchronized  with  the  same  or  a similar  time  transmission.  Preferably, 
for  greater  time-delay  resolution,  both  transmitting  and  receiving  sites  should  be  equipped 
with  atomic  time  or  frequency  standards.  When  the  transmitted  pulse  x(t)  is  of  short  dura- 
tion, the  sounder  receiver  essentially  measures  the  effective  linear  unit  impulse  response  func- 
tion h(t)  of  t'  at  channel.  The  received  signal  y(t)  is  given  by  the  convolution  integral: 

r °° 

y(t)=  I h(u)x(t  -uidu,  (30) 

— OO 


where  u is  a time  variable.  If  x(t)  is  an  approximate  impulse,  then  y(t)  is  obviously  propor- 
tional to  the  impulse  response  of  the  channel.  The  ettect  of  rapid  nuctuations  in  the  mode 
properties  can  be  reduced  if  several  pulses  are  transmitted  on  each  frequency  and  the  received 
responses  averaged.  Because  sounding  receivers  use  narrow  pulses,  their  bandwidths  are  com- 
monly 30  to  100  kHz;  therefore,  sounding  receivers  tend  to  be  more  susceptible  to  interfer- 
ence than  normal  hf  communications  equipment. 


r 


A most  important  variant  ol'  the  basic  pulse  sounding  system  is  one  in  which  the 
individual  pulses  on  each  frequency  are  modulated.  This  conlers  three  important 
advantages: 

1 . It  allows  pulse  compression  coding  to  be  applied  in  order  to  improve  the  time 
resolution  properties  of  the  system  without  having  to  shorten  the  basic  pulse 
length  and  hence  increase  the  peak  transmitter  power. 

2.  A limited  amount  of  data,  such  as  equipment  status,  can  be  transmitted  via  the 
pulse  modulation  coding. 

3.  It  improves  peuormance  in  high  interlerence  environments. 

Several  forms  of  pulse  compression  or  pseudorandom  codes  have  been  used  in  practice. 

These  include  Barker  codes,  Huffman  sequences,  and  binary  codes  (length  >13)  having 
approximately  aperiodic  autocorrelation  tunctions.  In  all  the  above  cases,  impulse  responses 
evaluation  is  accomplished  by  computing  the  input-output  cross-correlation  function  (ccf), 
0j^y(r),  as  follows: 

T/2 

f x(t)y(t+r)dt 
h(u)  - u)du  , 


(31 ) 

(32) 


where  rand  u are  time  delay  variables,  T is  the  correlation  interval,  and  is  the  input 

autocorrelation  function.  It  is  seen  from  equation  (32)  that  if  the  input  aperiodic  autocor- 
relation function  is  impulsive,  then proportional  to  the  system  impulse 
response  function. 

Coll  and  Storey' ' describe  an  application  of  pulse  compression  signals  to  the  field 
of  ionospheric  sounding.  Tliey  compare  single  pulses,  binary  codes  of  length  1 7,  and  a 
Huffman  impulse-equivalent  pulse  train  of  length  1 7.  The  basic  pulse  length  used  was  25  ps. 
The  energy  of  the  binary  code  relative  to  that  of  a single  component  pulse  with  the  same 
peak  power  was  1 7 (1 2.3  dB);  its  peak-to-maximum-side-lobe  ratio  was  8.5  to  1 (9.3  dB). 
The  Huffman  pulse  train  has  only  one  component  pulse  with  magnitude  of  1 and  hence  has 
less  energy  than  the  constant  magnitude  binary  signal.  The  energy  of  the  Huttman  signal 
was  7.28  (8.6  dB)  times  that  of  a single  pulse  of  the  same  peak  amplitude;  its  peak-to- 
maximum-side-lobe  ratio  was  10.8  (10.3  dB). 

Another  variation  of  the  pulse  sounder  is  a new  modulation  technique  based  on  the 

use  of  complementary  sequences.'  ^ if  one  of  the  pair  of  complementary  sequences,  say 

i 

I 

1 


* ' DC  Coll  and  JR  Storey,  Ionospheric  Sounding  Using  Coded  Pulse  Signal,  Radio  Science  Journal  of  Re- 
search, v 68D,  p 1155-1159 

I 

Darnell,  Channel  Estimation  Techniques,  paper  presented  at  NATO  AGARD  Conference  Proceedings: 
Radio  Systems  and  the  Ionosphere,  Athens,  Greece,  26-30  May  1975 
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X|(t),  is  applied  to  this  system,  it  will  give  rise  to  an  output  y](t).  Similarly  the  second 
sequence  x^(t)  will  yield  an  output  y-»(t).  Employing  (32)  for  each  signal  in  turn: 


I'l 


and 


. OO 


(33) 


(34) 


Summing  (33)  and  (34), 


,y,(T)=  J h(u)[^ 

“ “ OO  I- 


^x,x,‘^-u»+'^X2X2*^-'' 


du 


(35) 


l-quation  (35)  shows  that  the  sum  of  the  input-output  cross-correlations  at  equivalent 
values  of  t is  equal  to  the  convolution  of  the  system  impulse  response  and  the  sum  of  the 
two  individual  input  autocorrelation  functions.  This  assumes  that  the  two  sequences  are 
applied  to  the  linear  system  in  such  a way  that  there  are  no  cross-correlation  terms  between 
the  two  sequences. 

In  the  particular  case  of  ionospheric  sounding,  there  are  two  ways  in  which  the 
sequences  can  be  transmitted  in  order  to  fulfill  this  condition.  The  first,  illustrated  in 
figure  7a,  is  to  transmit  both  simultaneously,  in  time  synchronism,  as  near  to  each  other  in 
frequency  as  possible  without  their  spectra  overlapping.  The  respective  channel  cross- 
correlation responses  could  then  be  summed  to  give  the  desired  output.  A disadvantage  of 
this  mode  of  transmission  is  that  to  achieve  the  required  time  resolution  at  the  correlation 
output,  the  clock  rate  for  the  sequences  has  to  be  approximately  10  kHz  and  hence  the 
bandwidth  of  each  sequence  would  be  relatively  large.  Thus,  the  frequency  separation 
**  A ~ necessary  would  be  of  the  order  of  tens  of  kHz,  which  implies  that  propagation 
and  interference  conditions  would  differ  considerably  for  the  two  sequences.  In  the  second 
method,  illustrated  in  figure  7b,  the  complementary  sequences  can  be  transmitted  on  the 
same  frequency,  f^g,  sequentially,  with  sufficient  time  between  them  to  allow  for  the  decay 
of  all  multipath  returns  due  to  the  first  sequence  before  the  second  sequence  occurs  (typi- 
cally 10  ms).  This  again  means  that  there  is  no  cross-correlation  between  the  two  sequences. 
In  terms  of  processing  this  means  that  the  received  data  corresponding  to  the  first  sequence 
must  be  stored  and  correlated  in  synchronism  with  the  second  sequence.*  2 

Many  individual  pulse  sounders  are  in  operation  all  over  the  world.  The  US  Navy 
operates  the  Navy  Tactical  Sounder  System  at  three  transmit  sites.  France  operates  an 
identical  transmitter.  These  transmitters  transmit  two  13  bit  Barker  coded  sequences  25  ms 
apart  on  the  same  frequency  during  each  sweep.  F.ach  of  the  subpulses  is  of  200  ps  dura- 
tion. A 1 3 bit  Barker  coded  pulse  increases  the  received  signal-to-noise  ratio  by  13(11.1 
dB).  Hansen  and  Blitch*^  have  shown  that  the  Navy  Tactical  Sounder  receiver  has  a 


* '^Naval  Flfctronics  Laboratory  Center  TN  2399,  Modifications  to  the  AN/UPR-2  Sounder  Receiver,  by 
PM  Hansen  and  JD  Blitch,  13  June  1973.  Technical  Notes  are  working  documents  and  do  not  represent 
an  official  policy  statement  of  NOSC. 
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a.  PARALLEL 


FREQUENCY 


Figure  7.  Complementary  sequence  transmission 
and  processing.'  “ 

peak-to-maximum-side-lobc  ratio  of  2 (3dB)  without  the  use  of  an  interference  inhibitor 
circuit  and  5 (7.0  dB)  with  the  use  of  the  interference  inhibitor.  After  they  modified  tlie 
test  set  and  realigned  the  receiver,  they  obtained  a peak-to-maximum-side-lobe  ratio  of  1 2 
(10.8  dB).  They  designed  a new  pulse  amplifier  for  the  receiver;  its  peak-to-maximum-side- 
lobe  ratio  was  28  to  1 (14.4  dB). 

The  only  advantage  taken  of  the  two  sequences  in  the  NTS  receiver  is  in  averaging. 

If  a new  receiver  design  were  to  incorporate  a minicomputer  or  microprocessor,  it  would  be 
possible  to  store  the  first  arriving  Barker  coded  pulse  and  correlate  it  in  synchronism  with 
the  second  Barker  coded  pulse.  Figure  7b  illustrates  such  a process.  This  configuration 
would  provide  additional  processing  against  impulse  noise. 
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Figure  8 is  a tunctional  hloek  diagram  ot  a recommended  receiver  subsystem.  In 
addition  to  storing  the  first  transmitted  pulse,  the  microprocessor  could  be  programmed 
to  sense  what  is  happening  in  the  connected  receiver  and  could  issue  control  signals  in 
response  to  the  sensed  situation.  Hence  better  accuracy  in  data  collection  and  a fully  doc- 
umenteil  record  of  an  experiment  coukl  be  obtained.  To  insure  adequate  time  delay  resolu- 
tion. the  tnne-frei|uency  standards  should  be  atomic.  The  terminal  is  used  to  control  the 
processor  and  may  be  combined  with  the  optional  CRl  display  unit.  The  go/no-go  display 
IS  a cxmvenient  means  of  showing  ionospheric  conditions.  The  magnetic  tape  can  be  used 
as  a permanent  record  of  the  experiment.  Since  it  can  also  be  used  in  the  maintenance  ol 

the  equipment,  the  mag  tape  unit  should  be  both  a read  and  write  unit.  The  auxiliary  device  ! 

represents  oftline  equipment  not  under  direct  control  of  the  processing  unit.  This  device 

tiiay  be  a communication  system  needing  information  from  the  sounding  receiver,  it  may  be  | 

another  computer  processor,  or  it  may  be  a related  experiment,  l or  an  operating  wideband  | 

modem,  the  auxiliary  device  might  provide  the  necessary  frequency  dispersion  information  j 

neeiled  by  the  modem  for  its  proper  operation.  The  test  signal  emitted  by  the  sounding 

receiver  is  used  to  calibrate  the  receiver.  (That  is,  the  frequency  synthesizer  and  timing  of 

the  frequency  sweep  in  the  receiver  are  utilized  by  the  antenna  interface  unit  to  provide  the 

calibration.)  The  processor  controls  both  the  sounding  receiver  and  antenna  interface  unit. 

Another  form  of  ionospheric  sc.mding,  fundamentally  different  in  principle  from 
pulse  sounding,  is  "chirp"  sounding.  Again,  a synchronized  transmitter  and  receiver  are 
required;  now  . however,  a low  level  freciuency  sweep  is  transmitted  instead  of  pulses.  This 
sweep  can  typically  be  linear  or  logarithmic  in  form  and  can  cover  part  or  the  whole  ot 
the  hf  band.  ,\  replica  of  the  transmitted  sweep  is  generated  at  the  receiver  and  is  used 
effectively  as  the  local  oscillator  input.  The  signal  received  via  the  ionospheric  path  is  mixed 


f igure  8.  Rccoinmcndcd  receiver  subsystem  functional  block  diagram. 
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with  this  receiver  sweep,  anil  the  resulting  ilil't'erence  signal  is  subjecteil  to  spectrum  analysis. 
Multiple  paths  give  rise  to  multiple  receiveil  sweeps  with  various  amplitiules  and  relative 
delays.  Spectrum  analysis  sliows  up  tlie  dilTerence  I'reiiuencies  between  received  sweeps 
and  the  local  oscillator  sweep,  with  I'reiiuency  displacement  being  proportional  to  the  time 
delay  tor  a linear  sweep  and  amplitudes  being  correctly  translated  into  the  Ireiiuency 
domain.  A disadvantage  is  the  confusion  of  Doppler  with  time  delay. 


CLASS  2 CHANNEL  MEASUREMENTS 

The  channel  stimulus  .\(t)  is  transmitted  over  the  channel,  and  energy  which  is  re- 
llected  back  from  the  channel  in  tlie  form  of  a signal  y(t)  is  received  at  tlie  base  and 
used  in  the  channel  estimation  processes. 

Th  e processing  techniques  employed  in  practical  Class  II  systems  a^e.  for  tlie  most 
part,  very  similar  to  those  used  in  Class  I oblique  sounders;  again  they  essentially  evaluate 
a linear  impulse  response  function  for  the  channel.  However,  the  fundamental  distinction 
is  that  the  impulse  response  is  measured  solely  from  the  base  station.  Two  common  forms 
of  single-site  sounding  exist:  vertical  incidence  sounding  and  backscatter  sounding. 

A vertical  incidence  sounder  emits  a sounding  signal,  able  to  take  any  of  the  forms 
previously  described  for  oblique  sounding,  which  is  directed  vertically  at  the  ionosphere. 

The  retlected  returns  are  analy/.ed  by  a synchronized  sounding  receiver,  and  a vertical  incid- 
ence ionogram  is  produced.  Since  transmitter  and  receiver  are  now  collocated,  it  is  essential 
that  the  receiver  be  desensitized  during  transmission.  The  vertical  incidence  data  can  be 
used  directly  in  model  formulation  for  short  distance  skywave  paths;  for  long  range  paths, 
geometricai  computations  can  be  applied  to  the  ionogram  data  to  predict  oblique  incidence 
properties.  Wright  has  built  a prototype  digital  ionosonde  that  uses  an  associated  mini- 
computer and  which  has  very  high  time  delay  resolution.*'^ 

Backscatter  sounding  is  similar  to  vertical  incidence  sounding  except  that  the  trans- 
mitted energy  is  radiated  obliquely  instead  of  vertically.  The  received  signal  y(t)  is  now 
due  to  that  energy  which  has  propagated  via  the  ionosphere  over  a path  away  from  the 
base  and  has  been  reflected  from  the  ground  back  to  the  base  by  a similar  ionospheric 
reflection  in  the  reverse  sense.  Clearly,  attenuation  in  such  a system  is  large  and  much 
higher  transmitter  powers  are  required;  also  the  definition  in  the  backscatter  ionogram 
is  inferior  to  that  in  oblique  or  vertical  incidence  ionograms.  However,  backscatter 
sounding  does  give  a good  indication  of  usable  frequency  ranges  for  single-hop  paths; 
the  multihop  case  tends  to  introduce  excessive  path  loss. 


*^GH  Barry  and  RB  Fenwick,  Extraterrestrial  and  Ionospheric  Sounding  with  Synthesized  Frequency 
Sweeps,  Hewlett-Packard  J.  v 16.  p 8-12.  1965 

* ^National  Oceanic  and  Atmospheric  Administration  FRL  SHI  Preprint  20(i.  Development  ol  Sy,stems 
for  Remote  Sensing  of  Ionospheric  Structure  and  Dynamics:  Functional  Characteristics  and  Appli- 
cations of  the  “Dynasonde,"  by  JW  Wright.  July  1 975 
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CLASS  3 CHANNEL  MEASUREMENTS 

The  channel  stimulus  x(t)  and  the  base  data  input  bit)  are  multiplexed  prior  to 
transmission  over  the  hf  channel.  At  the  remote  station,  the  received  signal  y(t)  is  demul- 
tiplexed, processed,  and  used  to  t'ormulate  a model  which  then  determines  the  signal 
processing  strategy  to  be  applied  to  the  distorted  data  signal  b'it)  to  provide  the  estimate- 
bit). 

In  the  pilot  tone  technitiue,  low  level  cw  tones  are  either  inserted  in  the  data  stream 
or  transmitted  in  potentially  available  alternate  channels.  At  the  receiver,  simple  phase- 
measurements  on  the  tones  enable  the  relative  states  of  the  channels  to  be  specified  in 
terms  of  predicted  data  error  rates. 


TYPES  OF  SIMULATION  EXPERIMENTS 

There  are  several  approaches  in  using  a channel  simulator  to  measure  the  performance- 
capabilities  of  a modem.*  On  first  consideration,  the  usual  thought  is  to  select  a set  of 
typical  ionospheric  channels  representing  different  frequencies  in  the  hf  band  and  different 
transmitter-to-receiver  link  distances,  time  of  day,  latitude,  etc,  then  to  measure  the  modem 
performance  for  each  channel.  This  “typical  channel”  approach  will  show  the  modem’s 
capability  for  the  specific  combinations  of  channel  parameters  that  are  selected,  and  the 
results  can  be  compared  with  corresponding  measurements  on  the  other  modems.  It  is  in- 
efficient, however,  because  the  results  usually  cannot  be  used  to  determine  ( 1 ) the  contri- 
bution of  each  of  the  types  of  channel  distortion  to  the  modem’s  degraded  performance, 

(2 ) the  effects  of  modem-receiver  distortions,  (3 ) the  manner  in  which  the  various  distortions 
combine,  or  (4)  the  modem’s  performance  under  channel  conditions  for  which  measurements 
have  not  been  made. 

In  order  to  utilize  channel  simulation  to  its  finish,  it  is  necessary  to  apply  some 
significance  to  the  types  of  signal  distortions  caused  by  the  channel.  For  instance,  if  it  is 
found  that  a modem  under  study  has  a bit  error  rate  of  10~*  for  a Doppler  frequency  shift 
of  10  M/.,  it  is  necessary  to  know  how  and  under  what  circumstances  that  type  of  Doppler 
might  occur. 

A second  approach  that  could  overcome  these  inefficiencies  would  be  to  obtain  a 
complete  set  of  measurements.  In  the  “complete  set”  approach,  a suitable  number  of 
numerical  values  would  be  used  for  each  channel  parameter,  and  measurements  would  be 
made  for  all  combinations  of  parameter  values.  By  examining  selected  data  subsets  in  each 
of  which  only  one  distortion  parameter  changed  value,  the  inefficiencies  of  the  typical- 
channel  approach  could  be  overcome.  Unfortunately,  the  complete-set  approach  is  not 
practical  because  of  the  formidable  number  of  parameter  combinations. 

A compromise  approach  that  is  both  efficient  and  practical  might  be  called  the 
“single  and  multiple  distortions”  approach.  A channel  with  only  one  type  of  distortion  is 
used  first,  and  a set  of  measurements  of  bit  and  block  error  probabilities  are  measured  by 
using  successively  larger  values  of  the  distortion  parameter.  The  set  of  error  probability 
values  obtained  in  this  experiment  then  defines  bit  and  block  error  probability  curves  as 
functions  of  the  single  channel  distortion.  Similar  experiments  are  performed  for  other 
selected  channel  distortions.  The  experiments  are  then  expanded  to  include  two  types  of 
channel  distortion  at  a time,  and  error  probabilities  are  measured  as  functions  of  one 
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distortion  parameter  wliile  the  other  is  lield  constant.  The  procedure  is  continued  in  selected 
experiments  with  three  types  ot  cliannel  distortion  present,  where  error  probabilities  are 
always  measured  against  one  varied  channel  distortion  parameter  while  the  others  are  held 
constant.  Such  performance  curves  are  characteri/.ed  by  from  two  to  five  intervals.  With 
decreasing  abscissa  distortion  values,  the  curves  have  an  upper  plateau  interval,  followed  by 
at  least  one  descending  interval,  followed  in  most  cases  by  a second  plateau  interval,  and  in 
some  cases  by  a second  descending  interval  and  a third  plateau  interval.  In  each  interval,  one 
type  of  channel  or  modem-receiver  distortion  usually  predominates  and  the  other  distortions 
are  relatively  insignificant  in  determining  error  probabilities.  In  a small  fraction  ot  the 
intervals,  two  types  of  distortion  cun  be  comparable  and  significant.  Supplementary  measure- 
ments of  error  probabilities  versus  rms  signal  level  to  the  modem  receiver  can  be  made  with 
constant  channel  distortion! s),  to  yield  additional  information  about  modem-receiver 
distortions.  When  suitable  combinations  of  types  of  channel  distortion  are  chosen  with 
appropriate  parameter  values.the  various  curves  are  related,  and  the  approach  is  both  efiicient 
and  practical.  It  is  efficient  because  the  results  can  be  used  to  determine  the  extent  to  which 
each  type  of  channel  distortion  contributes  to  a modem’s  degraded  performance,  the  effects 
of  modem-receiver  distortions,  the  manner  in  which  the  various  distortions  combine,  and  a 
modem’s  performance  under  some  channel  conditions  for  which  measurements  have  not  been 
made.  The  approach  is  practical  because  a surprisingly  effective  set  of  measurements  can  be 
made  in  a reasonable  time.  It  is  the  method  that  was  used  in  the  channel-simulator  measure- 
ments on  the  USC-10,  ACQ-6,  and  MX-190  modems.'  The  specific  types  of  experiments  that 
were  performed  are  listed  in  table  2. 

TABLE  2.  TYPES  OF  CHANNEL-SIMULATOR  EXPERIMENTS.' 

Nonfading  experiments 
Signal  level 
Gaussian  noise 
Atmospheric  noise 
Doppler 

Fading  experiments 
Fading 

Signal  level  and  fading 

Gaussian  noise  and  fading 

Atmospheric  noise  and  fading 

Cw  interference  and  fading 

Doppler  and  fading 

Multipath  and  fading 

Fading  and  multipath 

Signal  level,  multipath,  and  fading 

Gaussian  noise,  multipath,  and  fading 

Atmospheric  noise,  multipath,  and  fading 

Cw  interference,  multipath,  and  fading 

Doppler,  Gaussian  noise,  and  fading 

Doppler,  multipath,  and  fading 
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An  hf  channel  simulator  provides  access  to  the  values  of  parameters  of  an  appropriate 
real-time  ionospheric  channel  model  which  describes  the  path  behavior  adequately. 

While  an  unproven  channel  model  can  provide  very  useful  information,  it  may  not  be 
typical  of  the  hf  medium.  A desirable  mathematically  tractable  channel  model  would  be  one 
that  is  known  to  be  a valid  representation  of  the  hf  medium  for  the  largest  practicable 
percentage  of  typical  hf  links. 

In  general,  there  are  two  basic  sources  of  signal  distortion  associated  with  hf  paths: 
time  dispersion  and  frequency  dispersion.  For  narrow  bandwidths,  the  ionospheric  delay 
variation  is  small  relative  to  the  multipath  spread,  and  dispersion  consequently  can  be  ignored 
in  simulating  hf  channels.  For  wide  bandwidths,  the  effect  of  multiple  modes  reduces  to  the 
production  of  reverberation.  Temporal  variations  of  the  ionosphere  no  longer  cause  fading, 
but  instead  produce  only  small  Doppler  shifts  (on  the  order  of  a fraction  of  a Hz)  between 
various  components  of  the  reverberation.  For  the  wide  bandwidths,  ionospheric  delay 
variation  can  cause  considerable  phase  distortion  in  which  the  phases  of  the  various  frequency 
components  of  a signal  are  altered  so  that  they  are  not  coherently  combined  by  the  receiver 
correlator. 

The  best  known  hf  channel  simulator,  developed  at  ITS  by  Watterson,  is  presently 
capable  of  accurately  simulating  all  types  of  channel  distortion  except  manmade  noise,  for 
system  bandwidths  less  than  2.5  kHz  at  night  and  less  than  1 2 kHz  during  the  day.  It  is  a 
four-tap  delay  line,  tach  tap  delivers  a delayed  but  effectively  undistorted  version  of  the 
transmitted  signal  to  a path  gain  multiplier.  Here  the  signal  is  modulated  in  amplitude  and 
phase  by  an  independent,  stationary,  complex,  baseband  function  termed  the  path  gain 
function,  Gj(t),  to  provide  the  frequency-scatter  distortion  (both  Doppler  and  fading)  on  that 
path. 

For  the  analysis  of  the  performance  capabilities  of  modems  suffering  multiplicative 
channel  distortions,  it  is  useful  to  have  statistical  descriptions  of  the  channel  frequency- 
scatter  and  time-scatter  characteristics  in  addition  to  the  channel-simulator  parameters. 
Watterson  derived  both  channel  frequency-scatter  and  time-scatter  characteristic  functions  in 
terms  of  the  path  gain  function  Gj(t).  The  derived  functions  include  channel  time-varying 
frequency  response,  H(f,  t);  channel  correlation  function,  R(Af,  At);  channel  scatter 
function,  S(t,  v)-.  channel  frequency  scatter  function,  V^.  (r*);  channel  attenuation,  a^,;  channel 
frequency  spread,  2a^.;  channel  time-scatter  function,  u^,  (r);  channel  delay,  and  channel 
time  spread,  2p^.. 

A source  of  distortion  of  ionospherically  propagated  signals  is  the  dispersion  caused 
by  ionospheric  refraction.  If  the  phase  distortion  imposed  on  a signal  is  represented  by  a 
transfer  function  of  the  form  exp(-i0(f)),  the  group  delay,  t,  becomes  d0/df.  If  r is  not 
constant,  then  dispersion,  with  its  consequent  distortion,  occurs.  The  dispersionjnay  be 
regarded  as  significant  if  b/2  (half  bandwidth)  is  equal  to  or  greater  than  ttIVAt/AL  where 
Ar/Af  in  s/Hz  is  the  slope  on  the  ionogram  at  f^,  the  center  frequency.  For  a 100  kHz 
bandwidth,  for  instance,  the  dispersion  will  be  significant  if  Ar/Af  is  equal  to  or  greater  than 
.TI4  X IO"^ms/kHz.  The  dispersion  tnay  be  represented  by  a filter  on  each  tap.  It  can  be 
approximated  by  a polynomial  function  of  frequency  at  a set  of  orthogonal  functions.  To 
measure  the  dispersion  accurately  enough  to  provide  the  coefficients,  the  ability  to  resolve 
time  delay  to  62.8  ps,  assuming  b = 100  kHz,  will  be  necessary.  The  time  frequency  scatter 
functions  must  be  recalculated  for  each  realization  of  time  delay. 


27 


i 


1 


! 


r 


In  Watterson  reported  an  experimental  verification  of  his  channel  model  in 
which  a hybrid  time-and-freiiuency  domain  method  was  used  to  obtain  propagation 
measurements.  The  signal  transmitted  was  amplitude  modulated  to  provide  coherent  pulses 
with  near-r.aussian  envelopes  of  50  /as  duration  at  200  llz.  Superimposed  on  this  pulse  were 
1 1 cw  tones  of  equal  amplitude.  I'he  1 1 tones  were  symmetrically  spaced  in  a kHz  band. 
Watterson  determined  that  his  model  was  valid  for  bandwidths  less  than  1 2 kHz  during  the 
daytime  and  2.5  kHz  during  the  nighttime.  The  sunrise  and  sunset  data  would  have  provided 
valuable  data  on  dispersion  if  the  dispersion  alone  had  been  analyzed.  Furthermore,  since  his 
transmitted  signal  consisted  of  a 50  /js  pulse  with  1 1 cw  tones  superimposed,  his  resolution 
for  time  delay  measurements  was  4.5  tis.  This  resolution  is  small  enough  for  the  necessary 
dispersion  correction  filters  to  be  determined. 

An  alternate  approach  to  hf  channel  simulation  given  here  is  useful  whenever  the 
ionospheric  dispersion  is  large  enough  to  cause  the  Watterson  model  to  require  an  excessive 
number  of  taps.  I'his  approach,  which  uses  fast  Fourier  transform  techniques,  provides  a 
transfer  function  for  a single  propagation  mode.  The  effects  of  multimoding  would  have  to 
be  handled  by  a superposition  of  such  transfer  timctions. 


RECO.Vl.MENDATIONS 


In  the  development  of  an  hf  channel  model  for  the  simulation  or  design  of  modems, 
include  the  effects  of  both  frequency  and  time  dispersion. 

To  obtain  representative  channel  model  data  necessary  to  implement  channel  simulation, 
use  an  oblique  sounder  system  that  achieves  needed  resolution  in  time  delay  by  employing 
complementary  sequences  or  modulated  pulses. 

If  the  Watterson  hf  channel  simulator  is  used  for  frequency  bandwidths  greater  than 
2.5  kHz  at  night  and  I 2 kHz  during  the  day.  represent  the  frequency  dispersion  by  a filter 
on  eacli  tap. 

If  the  spread  of  arrival  times  of  ionospherically  propagated  signals  would  cause  an 
excessive  number  of  taps  to  be  required  in  the  Watterson  model,  consider  the  alternate 
approach  of  using  fast  Fourier  transforms. 

Develop  the  receiver  subsystem  described  in  this  report  to  minimize  impulse  noise  by 
providing  additional  processing  gain. 

To  measure  the  performance  capabilities  of  a modem,  use  the  “single  and  multiple 
distortions"  approach  to  channel  simulation  described  in  this  report. 

To  determine  the  time  delay  versus  frequency  and  range  to  be  inserted  in  tlie  time 
delay  equation  for  each  mode  or  propagating  path,  use  the  Jones  ionospheric  ray  tracing 
program. 

(ienerate  a report  summarizing  the  channel  realization  intormation  reterenced 
herein  and  its  relationship  to  the  various  ionospheric  models  in  the  Jones  ray  tracing  program. 
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